Radiationless energy exchange in three-soliton collisions 
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We revisit the problem of the three-soliton collisions in the weakly perturbed sine-Gordon equation 
and develop an effective three-particle model allowing to explain many interesting features observed 
in numerical simulations of the soliton collisions. In particular, we explain why collisions between 
two kinks and one antikink are observed to be practically elastic or strongly inelastic depending 
on relative initial positions of the kinks. The fact that the three-soliton collisions become more 
elastic with an increase in the collision velocity also becomes clear in the framework of the three- 
particle model. The three-particle model does not involve internal modes of the kinks, but it 
gives a qualitative description to all the effects observed in the three-soliton collisions, including 
the fractal scattering and the existence of short-lived three-soliton bound states. The radiationless 
energy exchange between the colliding solitons in weakly perturbed integrable systems takes place 
in the vicinity of the separatrix multi-soliton solutions of the corresponding integrable equations, 
where even small perturbations can result in a considerable change in the collision outcome. This 
conclusion is illustrated through the use of the reduced three-particle model. 

PACS numbers: 



I. INTRODUCTION 

The study of soliton collisions in nonintegrable sys- 
tems 0, [1] is interesting because such systems typically 
describe more realistic situations than the integrable sys- 
tems where the interactions between solitons are known 
to be purely elastic [3J. In nonintegrable systems, the col- 
lision outcome can be highly nontrivial and, depending 
on the degree of nonintegrability, the collision scenario 
can have qualitatively different features. 



For the classical 4> equation, 
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which is rather far from an integrable system, kink col- 
lisions are always accompanied by a certain amount of 
radiation in the form of small-amplitude wave packets, 
as well as by the excitation of the kink's internal modes 
The latter are responsible for several effects in the 
(f> kink-antikink collisions. In particular, the resonant 
energy exchange between the translational motion of the 
kinks and their internal modes explains the fractal kink- 
antikink scattering [f|. This is a topic that was initiated 
by the numerical studies in Ref. @ (see also Ref. Q and 
references therein), and it is still under active investiga- 
tion 0. 

For long time, the excitation of the soliton internal 
modes and the radiation losses were thought to be two 
major manifestations of inelasticity of the soliton colli- 
sions in nonintegrable models. However, a qualitatively 
different manifestation was recently identified, namely, 
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the radiationless energy exch ang e (REE) between collid- 
ing solitons S[ISIi3Jp,IiililIiaIiEll3 (abbreviated 
as "REE" in Ref. jlOl] a designation that we will use 
hereafter) . 

The energy transferred to soliton internal modes in 
soliton collisions, for small e, is typically proportional 
to e 2 , and the same is true for the radiation losses (here 
e is the coefficient in front of a perturbation term, added 
to an integrable equation) . Terms proportional to e 2 ap- 
pear as the lowest-order correction terms in the collective 
variable approaches used to describe the soliton's inter- 
nal modes [J] ; the kink dynamics in the discrete <^> 4 equa- 
tion [lU ; the kink and breather d yna mics in the discrete 
sine-Gordon equation (SGE) jRj [2(| ; and the radiation 
from the discrete SGE kink [21| and from the soliton in 
the discrete nonlinear Schrodinger equation (NLSE) [1]. 
On the other hand, the degree of inelasticity due to the 
REE effect, when the latter is present (see details be- 
low) grows proportionally to e [lfl, EJ • This means that 
for weakly perturbed integrable systems the REE effect 
becomes a dominant manifestation of the inelasticity of 
collision, while the soliton's internal modes and radiation 
become increasingly important with stronger deviations 
from integrable case. 

The REE effect can also be responsible for the fractal 
soliton scattering which was demonstrated for the first 
time in [14| for the weakly perturbed SGE and later for 
the weakly perturbed NLSE [l5|, [l(| . In contrast to those 
studies, in Refs. [13, [HI the fractal scattering of vec- 
tor solitons in the coupled NLSE was attributed to the 
resonance energy exchange between the soliton's transla- 
tional and internal modes, i.e., through the mechanism 
similar to that operating for the 4 kinks [H, 0, 0, Q • 

Fractal soliton scattering in the weakly perturbed 
NLSE was explained qualitatively in the frame of a very 
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simple model [Hj] and for the generalized NLSE in the 
context of a more elaborate collective variable approach 
[2l | , based on the method of Karpman and Solov'ev [25[ . 
Remarkably, the soliton's internal modes were not in- 
volved into consideration in [lU, [24| indicating that the 
underlying dominant mechanism for the fractal scatter- 
ing was the REE effect (rather than the internal mode 
excitation) . 

For weakly perturbed integrable systems, parameters 
of the colliding solitons where the REE effect is observed 
can be found from the analysis of the corresponding inte- 
grable equation. This was done for the weakly perturbed 
SGE in [13] and for the weakly perturbed NLSE in ^ 
using the fact that the REE effect is observed in the vicin- 
ity of separatrix multi-soliton solutions of the integrable 
equation. 

In the case of moderate deviation from integrability, 
it becomes increasingly important to check if the degree 
of nonintegrability and the sign of perturbation allows 
for the existence of noticeable soliton internal modes be- 
fore one can judge to which extent the REE effect and 
the soliton internal modes contribute to the inelasticity 
of collision (see, e.g., Sec. II D in [l7[). The effect of 
the REE effect in the case of a moderate degree of non- 
integrability has studied far less extensively than in the 
case of weak perturbation, though valuable results have 
been recently obtained for the discrete NLSE [26] and for 
the generalized NLSE with various types of the nonlinear 
term [24| , where a general system of ordinary differential 
equations was derived for the velocities, amplitudes, po- 
sitions and phases of the solitary waves. The latter was 
shown to qualitatively and quantitatively match the pre- 
dictions of the full model. 

In this paper, we study the REE effect in three-soliton 
collisions of a weakly perturbed sine-Gordon equation. 
In the frame of the three-particle model, we demonstrate 
that the REE effect is directly related to a separatrix 
solution, and it offers a very transparent explanation of 
the origin of fractal soliton scattering. We also classify 
in a general way the potential for emergence of such phe- 
nomena in three-kink collisions of the weakly perturbed 
sine-Gordon model. 

The paper is organized as follows. In Sec. HH the 
results of numerical study of the degree of inelasticity of 
three-soliton collisions in the Frenkel-Kontorova model 
Eq. are presented. First, the collisions between three 
kinks/antikinks are analyzed in Sec. Ill Al and then the 
kink-breather collisions are investigated in Sec. Ill Bl The 
three-particle model is introduced and analyzed in Sec. 
IIIII The discussion of the results and our conclusions are 
presented in Sec. [TV] 



II. THREE-SOLITON COLLISIONS IN 
WEAKLY DISCRETE SGE 

To study the effects of non-integrability on the soliton 
collisions it is desirable to have a model with tunable 



deviation from an integrable case [c.f. with Eq. ([1} which 
does not have such a parameter] . The Frenkel-Kontorova 
(FK) model, 

-775 -^{Un-x - 2u„ + Un+i) + sinu„ = 0, (2) 

which is a discretization of the integrable SGE, 

iHt - u xx + sinu = 0, (3) 

is a convenient choice for such a study [13, EH [3 . The 
(singular) perturbation parameter in Eq. ^ is e = h 2 
(with h being the lattice spacing); the lowest order cor- 
rection to SGE due to the discretization can be quanti- 
fied, upon a Taylor expansion of the second difference, as 

The exact three-soliton solutions to SGE are well 
known [13, H3]. The solutions are the combinations of 
single-soliton solutions, namely kinks (K) or antikinks 
(K), having the topological charges q — 1 and q = — 1, 
respectively, and two-soliton solutions, namely breathers 
(B), which are actually the kink-antikink oscillatory 
bound states. 

Energy E and momentum P of one SGE kink are de- 
fined by its velocity V as follows 

E K = 85, P K = &V5, where tT 1 = y/l - V 2 . (4) 

Energy and momentum of a breather are defined by its 
frequency u> and velocity V: 

E B = Pb = 16^ 

where f _1 = y/l - V 2 , r) = yjl-ua 2 . (5) 

Below we describe the numerical results for the three- 
soliton collisions in the weakly discrete (h 2 = 0.04) SGE 
Eq. @. The exact three-soliton solutions to SGE were 
employed for setting the initial conditions. The equations 
of motion Eq. @ were integrated with the use of the 
Stormer method of order six. We register the parameters 
of quasi-particles after their collision and compare them 
with those before the collision. The larger the change in 
the parameters, the more inelastic the collision is. 

A. Three-kink collisions 

We number the kinks in a way that at t = (be- 
fore the collisions) their initial positions are related as 
(2:0)1 < (2:0)2 < (#0)3 and momenta as Pk x > 
Pk 3 - Here we consider only symmetric collisions with 
Pjc t > 0, Pk 2 — 0, and Pk 3 = —Pk x - Consideration 
of non-symmetric collisions does not bring any new im- 
portant physical effects. For the symmetric collisions it 
is convenient to set (2:0)1 = —(2:0)3 so that the three- 
soliton collisions are expected when (2:0)2 is close to the 
origin, otherwise the two successive two-soliton collisions 
will take place. Thus, among the kink's initial positions 
(2:0)1 the only essential parameter is (2:0)2- 
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Finally, our three-kink system is defined by the topo- 
logical charges of the kinks. There are eight possible 
variants in assigning the charges to the three kinks, 
which, due to symmetry, can be divided into three groups 
of topologically different collisions: KKK = KKK, 
KKK = KKK, and KKK = KKK = KKK = ~KKK. 
We will refer to each group by referring to their first mem- 
ber. 

The collision outcome is presented by the momenta of 
kinks after collision, Pr , as the functions of (2:0)2 f° r 
a given Pr 1 , which defines the initial momenta of the 
kinks, PXji a s described above. In some cases a kink- 
antikink pair can merge into a breather. In those cases we 
assumed that the kinks constituting the breather share 
its momentum equally, in order to plot their momenta. 

The results for the KKK collisions are shown in Fig. [T] 
(a) for Pr 1 — 0.8. Similar results for the KKK collisions 
are shown in Fig. [T](b) also for Pr 1 — 0.8. The results for 
the KKK collisions are shown in Fig. [2] (a) for Pr 1 = 2.5 
(larger collision velocity) and in Fig. [U(b) for P Kl = 0.8 
(smaller collision velocity). 

In the panels (a') and (b') of Fig. [1] and Fig. rj the 
examples of collisions are presented on the (x, t) plane 
by showing the regions of energy density greater than a 
certain value, so that the cores of the solitons are clearly 
seen. These examples are given for the particular val- 
ues of the initial coordinate of the middle kink, (#0)21 
indicated by the arrows in the corresponding panels at 
left. 

First we note that KKK and KKK collisions are al- 
ways practically elastic regardless of the specifics of (10)2 
(see Fig. [IJ and only KKK collisions are inelastic for 
(0:0)2 close to the origin (see Fig. [5]). We conclude that 
if a kink has positive or negative charge with equal prob- 
ability, then the REE in three-kink collisions can be ex- 
pected in two cases from eight. 

Of particular importance is the fact that in the strongly 
inelastic KKK collisions shown in Fig. [2] the energy 
given to the kink's internal modes and to the radiation is 
negligible in comparison to the energy exchange between 
the quasi-particles [l2j • This is the main feature of REE 
effect in soliton collisions. 

The KKK collisions can be strongly inelastic because 
in this case the cores of all three kinks can merge. Two- 
kink collisions are practically elastic for the considered 
case of weak perturbation, h 2 = 0.04, as it can be seen 
in Fig. [1] (b), (b'). To explain why the two-kink colli- 
sions are elastic we note that Eq. ([2]) conserves energy 
and, for small perturbation parameter h 2 , the momen- 
tum is also conserved with a high accuracy while the 
higher-order conservation laws of SGE are destroyed by 
the weak discreteness. The conservation of energy and 
momentum sets two constraints on the two parameters 
of the two-kink solution. A three-kink solution has one 
free parameter and REE becomes possible if all three 
kinks participate in a collision. 

For the KKK collisions we note that the collision with 
a larger velocity [see Fig. 2 (a), (a')] results only in quan- 
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FIG_1: Numerical results for (a), (a') KKK and (b), (b') 
KKK collisions in FK model. The left panels show the mo- 
menta of the kinks after collision Pk as the functions of the 
initial position of middle kink, (10)2- In both cases momenta 
of the kinks before the collision were Pk x = —Pk 3 = 0.8 and 
Pk 2 = and they are nearly same after the collision meaning 
that the collisions are practically elastic for any (3:0)2- The 
right panels show the examples of collisions on the (x,t) plane 
for (2:0)2 = by plotting the regions with the energy density 
greater than certain value, so that the cores of the solitons 
are clearly seen. 



titative change of kink parameters, while collision with a 
smaller velocity [see Fig. 2 (b), (b')] may result in fusion 
of a kink-antikink pair in a breather. The result of KKK 
collisions is extremely sensitive to variations in (2:0)2 m 
the vicinity of (2:0)2 = 0, especially for small collision 
velocities. 

A simple explanation of the fact that the collisions be- 
tween two kinks and an antikink are always practically 
elastic for KKK and can be strongly inelastic in the case 
of KKK will be offered in Sec. ifll 

In the case of weak perturbation we never observed 
fractal patterns in the three-kink collisions (recall that 
in the tfi 4 model such patterns can be observed even in 
two-kink collisions but, as it was already mentioned, this 
model is far from an integrable one), while it can be ob- 
served in the kink-breather collisions, as discussed below, 
and in the breather-breather collisions UM. 



B. Kink-breather collisions 

Without loss of generality, we assume Pk + Pb = 0. 
Then we have two parameters, the momentum Pg and 
frequency u of the breather. The outcome of the KB 



4 




FIG. 2: Numerical results for KKK collisions with (a), (a') 
Pk\ = 2.5 (larger collision velocity) and (b,b') P«,\ = 0.8 
(smaller collision velocity) in the FK model. The left panels 
show the momenta of the kinks after collision Pk as the func- 
tions of the initial position of middle kink, (20)2- Collisions 
are strongly inelastic for (10)2 close to the origin. The right 
panels show the examples of collisions on the (x, t) plane for 
(a') (£0)2 = 0.05 and (b') (2:0)2 = 0.01 by plotting the re- 
gions with the energy density greater than a certain value, 
so that the cores of the solitons are clearly seen. Collision 
with a larger velocity in (a), (a') results only in quantitative 
change of kink parameters while collision with a smaller ve- 
locity in (b), (b') may result in fusion of a kink-antikink pair 
in a breather. 



collisions is studied as a function of the initial separa- 
tion between the kink and the breather controlled by the 
initial kink position (xq)k ■ 

In Fig. [3] we plot the momenta of kinks after collision, 
Pkj (including the kinks constituting the breather, as- 
suming as earlier that they share the breather's momen- 
tum equally), as a function of (xq)k for (a) Pg = 2.5 
(larger collision velocity) and (b) Pb = 1-6 (smaller col- 
lision velocity). One can see that strong REE is possible 
in the KB collisions. Note that in Fig. [3] only a small 
part of one period of the output functions is shown for 
the region with strong REE. In (a) there is a range of 
(xq)k where the breather obtains enough energy to split 
into a kink-antikink pair [example is shown in (a')]. In 
(b), in addition to this possibility, there appears a region 
where the breather is reflected from the kink [example is 
shown in (b')]< 

For large collision velocities [somewhat larger than in 
Fig. [3] (a), (a')] the kink passes through the breather 
with no qualitative change in the collision outcome; there 
is only some energy and momentum exchange between 



FIG. 3: Numerical results for the kink-breather collisions in 
FK model for (a), (a') Pb = 2.5 (larger collision velocity) and 
(b,b') Pb = 1.6 (smaller collision velocity) and u = 0.05 in 
both cases. The left panels show the momenta of the kinks 
after the collision Pk ■ as a function of the initial position of 
the kink, (xq)k- The right panels show examples of collisions 
on the (x,t) plane for (a') (xo)k = 40.25 and (b') (xq)k = 
28.1. 



them. However, for sufficiently small collision velocity 
the collision outcome as the function of the initial sep- 
aration between the kink and breather is a fractal. An 
example is presented in Fig. 0] for uj — 0.3, Pb = 
(kink and breather have zero initial velocities), where 
the soliton's momenta after collision Pk are shown as 
the functions of (xq)k- 

The structure presented in Fig. [4] can be described as 
a chain of self-similar patterns. At each scale smooth 
regions are separated by the apparently chaotic regions 
of two symmetry types, one shown in (a) and (d) and 
another one in (b) and (c). (b) and (d) present blowups 
of the regions indicated in panel (a) ; (c) is a blowup of 
the region indicated in (b). 

Two examples of the kink-breather dynamics are given 
in Fig. [5] for (a) (a*)* = 0.24 and (b) (x ) K = 0.236 
[indicated in Fig. [1] (b) by the arrows A and B, respec- 
tively] . The three-particle solution has a certain lifetime 
L (in this example L « 55) and then it splits into a kink 
and a breather. Similar dynamics has been reported, e.g., 
for the breather-breather system in the weakly discrete 
FK model [3, in th e weakly perturbed NLSE and 
recently for the generalized NLSE [24]. Thus, this type of 
dynamics is rather general. For the two-soliton collisions 
in the weakly perturbed NLSE we have estimated nu- 
merically the probability p to observe the three-particle 
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FIG. 4: Fractal kink-breather scattering observed for a; = 
0.3, Pb = (kink and breather have zero initial velocities). 
The kinks' momenta after the collision Pk- are shown as the 

3 

functions of the initial position of the kink, (xo)k, at different 
scales. At each scale smooth regions are separated by the 
apparently chaotic regions of two symmetry types, one shown 
in (a) and (d) and another one in (b) and (c). (b) and (d) 
present blowups of the regions indicated in panel (a); (c) is a 
blowup of the region indicated in (b). 
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FIG. 6: Probability p to observe the kink-breather system 
with the lifetime L (in Fig. [5] we have L ~ 55). Numerical 
data is shown by dots and only the cases with L > 1000 were 
taken into account. Dashed line is the guide for an eye and it 
has slope -3.5. 



including the existence of the two qualitatively different 
stochastic regions in the fractal structure will be clarified 
in Sec. IIHI with the help of the three-particle model. 



III. THREE-PARTICLE MODEL 




FIG. 5: Examples of the kink-breather dynamics for (a) 
(x ) K = 0.24 and (b) {x )k = 0.236 [indicated in Fig. g] 
(b) by the arrows A and B, respectively]. 



system with the lifetime L and found that p ~ L~ 3 
[15j | . Here we carry out a similar estimation for the kink- 
breather solution in the FK model and the result is shown 
in Fig. [5] The numerical data can be fitted as p ~ L~ 3 5 . 
There is evidence that for sufficiently small frequency of 
the breather the kink-breather system in the FK model 
with a small h 2 never splits [l^ . 

All the important features of the KB fractal scattering 



A. Description of the model 

Attempting to explain the effects observed in the three- 
soliton collisions in weakly perturbed SGE reported in 
Sec. [HI we consider the solitary waves as effective par- 
ticles, and study the dynamics of three such particles in 
one-dimensional space. The particles have mass m = 8, 
which is the rest mass of SGE kink, and they carry topo- 
logical charges qj = ±1. Particles with qj = 1 (qj = — 1) 
will be called kinks (respectively, antikinks) by analogy 
with the SGE solitons. We assume that particles i and j, 
having coordinates Xi and Xj , interact via the potential 

Uij (Tij ) = 16 + j , m = - x, , (6) 

which qualitatively approximates the interaction of two 
SGE kinks. The potential of Eq. © is attractive for 
q-i ^ qj and repulsive for q^ = qj. The binding energy of 
the kink-antikink pair is equal to 16, which is the energy 
of two standing SGE kinks. Note that for the kink and 
antikink at any finite distance ry the potential energy 
Uij {rij) is less than 16. If the kinetic energy of relative 
motion of the particles is less than 16 — Uij (rij), then the 
particles cannot escape the mutual attraction and they 
form an oscillatory bound state, i.e., a breather. 
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The Hamiltonian of the three-particle system is 



H = 7 + ^12^12) + U 13 (ri 3 ) + C/23(r 23 ) , 

j'=i 



(7) 



where Uj = dxj/dt, and there is one more integral of 
motion, namely the conservation of momentum. Without 
loss of generality, we assume that the total momentum 
in the system is equal to zero, i.e., m(«x + v 2 + w 3 ) = 0. 
Introducing new variables 

X2 — x\ — > v3a + /3, x 3 — xi — ■> 2/3, i — > V2mt , (8) 

the Hamiltonian of Eq. ([7]) can be presented in the form 

H = 1 (d 2 + /3 2 ) + lq 2 (V3a + /3) 



+U 13 (20) + U 23 (V3a - 0) . 



(9) 



which is the Hamiltonian of a unit-mass particle moving 
in the two-dimensional scattering potential. 

Now we solve numerically three equations of motion 
which can be derived from the Hamiltonian Eq. {7} and, 
inverting Eq. ((SJ) , present the three-particle dynamics by 
the trajectory of the particle in the (a, /3)-plane. 
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FIG. 7: (Color online) (a) The SGE solution of Eq. (JTOJ for 
V — 0.2; (b) the three-particle dynamics with qi = —52 = 
93 = 1, («o)i = -(^0)3, («o)a = 0, (aso)i = -(3:0)3 = -25, 
and (0:0)2 = in the (3;, t) space; (c) the red line shows the cor- 
responding trajectory of the particle in the scattering poten- 
tial in the (a, /3)-plane (isopotential lines are shown in black). 
The particle in (c) moves along the potential ridge and this 
motion is unstable. The picture in (c) gives a visual image of 
the separatrix KKK solution Eq. (|10|l . 



B. Separatrix three-soliton solutions to SGE 

Several separatrix three-soliton solutions to the exactly 
integrable SGE Eq. ([3]) have been reported in [12| . Here 
we reproduce two solutions important for our study. 

The separatrix KKK solution is 

R 

u k~kk ( x 't) = ^ arctan(exp x) + 4 arctan — , 
R = 5 (sinh F — cosh G sinh x) , 
S — 8 (cosh G + sinh F sinh x) — cosh F cosh a;, 

F = -Sx, G = SVt, 5- 1 = y/l - V 2 . (10) 

In this highly symmetric solution the anti-kink is at rest 
and it is located at the point of collision of two kinks 
moving with the velocities V and — V. 
The kink-breather separatrix solution is 

u kb ( x Tt) =4 arctan(exp x) + 4 arctan — , 

X = rj (sinh D — cos C sinh x) , 
Y = f] (cos C + sinh D sinh x) — cosh D cosh a;, 

C = -Lot, D = i]x, r)= y/l -oj 2 , (11) 

and it has only one parameter us because it is a particular 
form of the KB solution where the kink and the breather 
have zero velocities and zero distance between them. 

In Fig. Ewe plot (a) the SGE solution Eq. (pU for V = 
0.2, (b) the three-particle dynamics in the (x, t) space for 

qi = -<72 = qs = l, {vq)i = -(u ) 3 , (v ) 2 = 0, (x )i = 




xx a 



FIG. 8: (Color online) (a) The SGE solution Eq. (JTTJ for 
us = 0.2; (b) the three-particle dynamics with qi = —52 = 
93 = 1, («o)i = (1*0)2 = (uo)s = 0, (3:0)1 = -(3:0)3 = -4, and 
(3:0)2 = in the (x,t) space; (c) the red line shows the corre- 
sponding trajectory of the particle in the scattering potential 
in the (a, /3)-plane (isopotential lines are shown in black). The 
particle in (c) oscillates along the potential ridge and this mo- 
tion is unstable. The picture in (c) gives a visual image of the 
separatrix KB solution Eq. 



— (3:0)3 = —25, and (2:0)2 = 0, and in (c) the red line 
shows the corresponding dynamics in the (a, /3)-plane. In 
(c) the isolines of the scattering potential are also shown 
(black lines). The scattering potential in this case is a 
superposition of a ridge along [3 = and two troughs 
along the lines [3 = ±v3cv. Note that the intersection of 
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the ridge and the two troughs forms in the vicinity of the 
origin the ridge along the line a = 0; the trajectory of 
the particle shown by the red line goes exactly on the top 
of this ridge. Obviously, this type of motion is unstable 
and, as we will see in the following, small variation in 
the initial conditions may result in qualitatively different 
dynamics of the particle. The picture presented in Fig. [7] 
(c) gives a visual image of the separatrix KKK solution 

Eq. unai. 

In Fig. ED we plot (a) the SGE solution Eq. {TTJ) for 
u) = 0.2, (b) the three-particle dynamics in the (x, t) 
space for q x = -q 2 = q 3 = 1, (v )i = (uo)a = (^0)3 = 0, 
( x o)i = — (^0)3 = — 4, and (^0)2 = 0, and in (c) the 
red line shows the corresponding trajectory in the (cy, /?)- 
plane. The particle in (c) oscillates along the potential 
ridge and, similarly to the previous example, this motion 
is unstable. The picture presented in Fig. [8] (c) gives a 
visual image of the separatrix KB solution of Eq. (jTTJ) . 

When the red line passes the origin of the (a, /3)-plane, 
from Eq. |8]) one has x\ — X2 = £3, i.e., all three parti- 
cles meet at one point. In the SGE this corresponds to 
simultaneous collision of all three kinks. 

Looking at Fig. [UJc) one can expect the possibility of 
oscillation of the particle along the ridge of the scattering 
potential with j3 = 0. This is indeed possible for the 
three-particle system but, from (3 = one finds from Eq. 
[5] that xi(t) = X3(t), which cannot be realized in the FK 
model because the kinks have finite width. 






KKK 




< 



96 P 




FIG. 9: (Color online) Comparison of (a,a') KKK and (b,b') 
KKK symmetric collisions. The top panels show the three- 
particle dynamics in the (x, t) space. The trajectories of kinks 
are shown by thicker lines than those of antikinks. The bot- 
tom panels correspondingly show the equipotential lines of 
the scattering potential Eq. © (black) and the trajectory 
of particle (red line) in the (a, /3)-plane. The parameters are 
(aco)i = -(2:0)3 = -25, (0:0)2 = and (u )i = -(^0)3 = 0.6, 
(1*0)2 = 0. The charges of the particles are (a,a') q\ = q% = 
qz = l and (b,b') gi = q 2 = -93 = 1- 



C. Three-kink collisions 

In Fig. [9] we compare the KKK and KKK sym- 
metric collisions in the three-particle model for (xq)i — 
-(2:0)3 = -25, (£0)2 = and (w )i = -(^0)3 = 0.6, 
(■^0)2 = 0- The top panels show the three-particle dy- 
namics in the (a;, t) space. For each case, the bottom 
panels correspondingly show the equipotential lines of the 
scattering potential of Eq. (black) and the trajectory 
of the particle (red line) in the (a, /3)-plane. In (a') the 
scattering potential for qi = q2 = q 3 = 1 is a superposi- 
tion of three troughs while in (b') for q\ = q2 = —93 = 1 
it is a superposition of a ridge and two troughs. The po- 
tential in Fig. [5] (b') can be obtained from that shown in 
Fig. [7| (c) and Fig. [H] (c) through a rotation by — 7r/3. 

In Fig. [U (a), the like particles repel each other and, 
in (a'), the particle hits the potential barrier and goes 
back. In (b), one can see that particles collide in two 
successive two-soliton collisions. In this case, the particle 
in (b') passes the two potential troughs one after another 
and then moves away from the origin in the direction 
symmetrically equivalent to the direction it came from. 
Since the red line in (a') and (b') never goes through the 
origin, the three particles never meet at one point. 

In Fig. [10] we give two examples of near-separatrix 
KKK symmetric collisions in the three-particle model 
for (x )i = -(afo)3 = -25 and (v )i = -(v ) 3 = 0.6, 
(^0)2 = 0- Recall that the separatrix solution shown in 
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FIG. 10: (Color online) Sensitivity of the result of near sep- 
aratrix collision to a small deviation from (0:0)2 = demon- 
strated by setting (a, a') (0:0)2 = 1.2 and (b,b') (0:0)2 = 0.2. 
In (a,a') only a quantitative change in the system can be seen 
upon collision [compare with the actual three-kink collision in 
FK model shown in Fig. H]V)]. In (b,b'), kink and antikink 
merge in a breather [compare with Fig. [2jb')]. Other pa- 
rameters: (a;o)i = —(0:0)3 = —25 and («o)i = —(^0)3 = 0.6, 
M2 = 0. 
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Fig. [7] corresponds to (2:0)2 = but Fig. [TOl corresponds 
to (a,a') (x ) 2 = 1.2 and (b,b') (x ) 2 = 0.2. In Fig. [TO] 
(a, a'), the deviation from the separatrix is rather large 
and only quantitative changes in the particle parameters 
can be seen. This should be compared with the actual 
three-kink collision in FK model shown in Fig. [2fa'). In 
(b,b'), kink and antikink merge in a breather [compare 
with Fig. |2fb')]. Taking into account the time reversibil- 
ity in the Hamiltonian systems, this picture can be also 
regarded as an illustration of the breakup of a breather 
colliding with a kink. 

The three-particle model explains why the REE effect 
is more pronounced for the solitons colliding with a small 
relative velocity. For the particle moving in the (a, /?)- 
plane along the separatrix [red line in Fig. 0c)], any per- 
turbation results in exponential in time deviation from 
the potential ridge. High-speed collision results in faster 
passing of the scattering potential and the trajectory of 
the particle cannot be considerably changed. The situa- 
tion is opposite for the slow particle, which corresponds 
to the collision of solitons with a small relative velocity. 

D. Kink-breather collisions 

Here we select the parameters of the three particles 
so as to simulate the collisions between a kink and a 
breather. In particular, we set the charges of particles 
as qi = —qi =#3 = 1, their initial velocities as (i>o)i — 
(^0)2 = 0.3, (1*0)3 = —0.6; the initial positions of the 
particles constituting the "breather" are (20)1 = —16, 
(^0)2 = —13.5, and the third particle initial position was 
varied. In Fig. [11] the results are shown for (a) (2:0)3 = 
30.51, (b) (2:0)3 = 23.398, and (c) (#0)3 = 23.391. The 
top panels show the three-particle dynamics in the (x, t) 
space, while the bottom panels show the corresponding 
trajectory of the particle in the (a, /3)-plane (red line). 

Collisions in (a) and (b) are elastic but the difference 
is that while in (a') the particle does not move along the 
separatrix line a = 0, in (b') it does, and a very small 
change in the initial conditions is sufficient to have a qual- 
itatively different result of the collision, as presented in 
(c),(c'), where the breather reflects from the kink [com- 
pare (c) with actual kink-breather collision in the FK 
model shown in Fig. 3(b')]. 

E. Fractal kink-breather scattering 

To reproduce the kink-breather fractal scattering de- 
scribed in Sec. Ill Bl for the FK model we set the fol- 
lowing parameters for the particles in the three-particle 
model: qi = -q 2 = q 3 = 1, (u )i = (1:0)2 = (^0)3 = 0, 
(2:0)1 = -(2:3)2 = -5, and variable (2:0)2- 

In Fig. [T2] we present the velocities of particles after 
collision Vj as the function of (2:0)2- In (b) a blowup of 
the self-similar region indicated in (a) is presented. Com- 
parison of the panels (a) and (b) with the corresponding 




FIG. 11: (Color online) Three-particle model simulating the 
kink-breather collisions. The top panels show the three- 
particle dynamics in the (x, t) space, while the bottom pan- 
els show the corresponding trajectory of the particle in the 
(a,/3)-plane (red line). Only the initial position of the third 
particle is varied: (a) (3:0)3 = 30.51, (b) (2:0)3 = 23.398, and 
(c) (0:0)3 = 23.391. Collisions in (a), (a') and (b),(b') are elas- 
tic but in the latter case it is close to the separatrix (see Fig. 
[5} resulting in a great sensitivity to variations in initial condi- 
tions, as demonstrated in (c),(c'). The rest of the parameters 
are chosen as gi = —92 = qa = 1, (i>o)i = (^0)2 = 0.3, 
(«o)3 = -0.6, (2-0)1 = -16, (20)2 = -13.5. 




0.41 0.42 , 0.43 0.44 



FIG. 12: Fractal three-particle scattering. Panels (a) and (b) 
should be compared with the corresponding panels of Fig. [4] 
Parameters: qi = -qi = q 3 = 1, (v )i = (1*0)2 = ("0)3 = 0, 
(2:0)1 = —(2:3)2 — —5, and variable (20)2- 

panels of Fig. [5] reveals the qualitative similarity in the 
KB collision outcome in the FK model and in the three- 
particle model. We note that while we expect this parti- 
cle model to bear the essential qualitative characteristics 
of the three-particle collisions, their details depend sen- 
sitively on the precise initial conditions; for this reason, 
we expect Fig. [T^] to match qualitatively the results of 

Fig. m 

One of the important elements of the usefulness of the 
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FIG. 13: (Color online) (a) The exact KB solution to the 
SGE, (b) the same solution in the weakly discrete FK model, 
(c) the three-particle dynamics in the (a;, i) space and (d) the 
corresponding dynamics in the (a, (3) space. In (a), (b), kink 
and breather have zero initial velocities, breather frequency is 
io — 0.3, and separation between the kink and the breather is 
equal to 1.2. In (c), (d), gi = -q 2 = q 3 = 1, («o)i = (^0)2 = 
(«o)3 = 0, (x )i = -6.18, (20)2 = 0, (20)3 = 3.3. 



three-particle model is that it gives the possibility to an- 
alyze the KB fractal scattering from a different point of 
view, namely, by looking at the corresponding dynamics 
of the particle in the scattering potential in the (a, (3) 
space. In Fig. [T3] we show (a) the exact KB solution 
to the SGE, (b) the same solution in the weakly discrete 
FK model, (c) the three-particle dynamics in the (x, t) 
space and (d) the corresponding dynamics in the (a, (3) 
space. In (a), (b), kink and breather have zero initial 
velocities, breather frequency is u> = 0.3, and separation 
between the kink and the breather is equal to 1.2 (we 
refer to the form of the KB solution given in [l2j ). In 
(c), (d), q 1 = -q 2 = q 3 = 1, («o)i = («o)2 = Oo)3 = 0, 
(2:0)1 = —6.18, (2:0)2 = 0, (2:0)3 = 3.3. One can see 
from Fig. [T3]that the distance between the kink and the 
breather does not change in time in the integrable sys- 
tem [shown in (a)] but in the nonintegrable ones the dis- 
tance between them gradually decreases and they even- 
tually collide [see (b) and (c)]. The separated kink and 
breather having zero velocities are presented in the (a, (3) 
space by the particle oscillating along the line normal to 
the trough with orientation (3 — \/3a (kink is to the left 
of the breather) or /3 = — Via (kink is to the right of 
the breather) . However, the troughs have a slope toward 
the origin of the (a, (3) plane and the oscillating particle 
gradually approaches the origin, i.e., the collision point 
of three particles. After the particle has approached the 
origin [this situation is shown in Fig. [131(d)], two qualita- 
tively different scenarios giving different fractal patterns 
are possible. 

The first scenario is shown in Fig. 1141 Here the parti- 
cle after making a few oscillations normally to the trough 
(3 = \/3a can cross the separatrix line a — and make 
a few oscillations normally to the trough (3 = — \/3a and 
then again cross the separatrix line changing the trough. 
While the particle performs such crossings of the sepa- 
ratrix it remains close to the origin of the (a, (3) plane 




a 



FIG. 14: (Color online) Illustration of one of the two possible 
scenarios of fractal kink-breather scattering. (a),(b) show the 
dynamics of three particles in the (x,t) plane while (a'),(b') 
the corresponding dynamics in the (a, fi) plane. The parti- 
cle in (a'),(b') after making a few oscillations normally to the 
trough P = \> r 3a can cross the separatrix line a = and make 
a few oscillations normally to the trough (3 = — v3a and then 
again cross the separatrix line changing the trough. While 
the particle performs such crossings of the separatrix it re- 
mains close to the origin and thus, all three particles are close 
to each other. This defines the lifetime of the multi-soliton 
bound state discussed in [3, [TBI |24| . Eventually, the parti- 
cle will move away from the origin along one of the troughs. 
Parameters: (vo)i = 0, i = 1,2,3, (xo)i = — (20)3 = —5, (a) 
(x ) 2 = 0.433 and (b) (350)2 = 0.4275. 



and thus, all three particles are close to each other. This 
defines the lifetime of the multi-soliton bound state dis- 
cussed in [3, EH, HH ■ The probability to have a three- 
particle bound state with a long lifetime is small (sec Fig. 
[6]) meaning that eventually the particle will move away 
from the origin along one of the troughs remaining in the 
half- plane (3 > (compare Fig. [14] with Fig. [5] where the 
KB dynamics in the FK model is presented). 

The second scenario is more obvious because it is di- 
rectly related to the separatrix solution Eq. (fTTj) pre- 
sented in Fig. [H After making a few oscillations nor- 
mally to the trough j3 — \/3a as shown in Fig. [T3] (d) , 
the particle can be sent by the scattering potential al- 
most exactly along the separatrix line a = 0. Then the 
particle will make several oscillations along the ridge of 
the potential, as shown in Fig. [5] (c), before the inherent 
instability of this trajectory "ejects" the particle away 
from the origin in one of the four directions along the 
troughs (3 — ±-\/3a. This contrast to the first scenario 
where the particle can be scattered by the potential in 
the two of the four directions, namely, in the ones with 
(3>0. 

The first scenario is associated with the self-similar re- 
gions connecting the two "butterflies" [see Fig. 2] (b) 
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and (c)] while the second one is associated with the self- 
similar regions connecting the "wings" of a "butterfly" 
[see Fig. [4] (d)]. However, the whole fractal pattern is 
the result of the combination of both mechanisms. Each 
scenario is related to a periodic orbit of the particle in the 
scattering potential [28J . In the second scenario the peri- 
odic orbit is the separatrix kink-breather solution shown 
in Fig. [51 while in the first scenario there exists an in- 
finite set of periodic orbits. One particular orbit can be 
described as follows: the particle in the scattering poten- 
tial makes N oscillations almost normally to the trough 
ft = ^/3a and then jumps to the trough ft — —y/3a 
where it also makes N oscillations and then returns to the 
trough ft = y/lia completing one period of the periodic or- 
bit [Fig. [TJ] (a') and (b') give examples when the particle 
makes such jumps between the troughs ft — ±V3a but in 
these cases the trajectories are aperiodic] . If the particle 
follows a periodic orbit exactly, the three-particle sys- 
tem never experiences a breakup; however, the eventual 
separation of the structures is a result of the dynamical 
instability of such periodic orbits. 

It is well-known that the probability p of the time delay 
T for the particle interacting with the scattering poten- 
tial without the periodic orbits decreases exponentially 
with T while in the presence of the periodic orbits it de- 
creases algebraically [28| . The scattering potential in our 
case does have the periodic orbits and the probability p 
to observe a bound state with the lifetime L (analogous 
to the time delay T) decreases algebraically, p ~ L~ a . 
This was found in [15] for the two-soliton collisions in 
the weakly perturbed NLSE, and in the present study 
this was also confirmed for the kink-breather system in 
the FK model, as presented in Fig. [51 



IV. DISCUSSION AND CONCLUSIONS 

Through direct numerical simulations, we have pre- 
sented some of the striking effects generated by even a 
weak breaking of integrability (via discretization) in the 
sine- Gordon model. We have indicated that alternative 
mechanisms such as the excitation of internal modes and 
the emission of phonon radiation are too weak to ex- 
plain the phenomena observed in numerical simulations, 
and we have therefore attributed them to the radiation- 
less energy exchange between the solitons. Indeed, these 
effects have been systematically explained in a qualita- 
tive fashion in the framework of the three-particle model 
suggested in Sec. IIII Al lending direct support to the 
conclusion that all the nontrivial effects are due to the 
radiationless energy exchange between colliding solitons 
@, ©, M El HI M M HE HE E3- The following is 
known about the REE effects: (i) Manifestations of the 
REE effect grow proportionally to the perturbation pa- 
rameter e while radiation and excitation of soliton inter- 
nal modes grow as e 2 . (ii) In the sine-Gordon model the 
REE effect can happen only if at least three solitons col- 
lide simultaneously. Energy exchange in the two-soliton 



collision is suppressed by the two conservation laws that 
remain exactly or approximately preserved in the weakly 
perturbed system, (iii) The REE effect is related to 
the existence of the separatrix multi-soliton solutions to 
the integrable equations. Near-separatrix motion is ex- 
tremely sensitive to the perturbations [29] . (iv) The REE 
effect can be responsible for the fractal soliton scattering. 

The REE effect is generic and some of the above con- 
clusions can be also extended to other nearly integrable 
models 0]. For instance, the REE effect is observed in 
the weakly perturbed NLSE already in two-soliton colli- 
sions because here each soliton has two parameters and 
the total number of parameters describing the two-soliton 
solution (four) exceeds the number of the remaining con- 
servation laws. On the other hand, the REE is not pos- 
sible in the weakly perturbed KdV equation or weakly 
perturbed Toda lattice [30j because in these cases the 
soliton's cores never merge during collisions and thus, 
the multi-particle effects arc absent. Interestingly, the 
fractal pattern of different nature (not related to REE) 
is possible in KdV systems [31 1. 

The three-particle model offered in the present study 
can be reduced to the study of the dynamics of a par- 
ticle interacting with the two-dimensional scattering po- 
tential. Such a reduction gives a clear interpretation of 
the abovementioned features of the REE effect observed 
in the three-soliton collisions. Particularly, the following 
features have been identified: 

1. The three-particle model gives a visual image of the 
separatrix three-kink and kink-breather solutions 
to the integrable SGE, see Fig. [7J and Fig. [SJ The 
separatrix solution corresponds to motion of the 
particle along a ridge of the scattering potential. 

2. KKK collisions and KKK collisions are al- 
ways practically elastic while KKK collisions are 
strongly inelastic in the vicinity of (xo)2- Only in 
the latter case the point moves along the ridge of 
the scattering potential, which is the motion along 
a separatrix, see Fig. [Jj Also only in KKK colli- 
sions the point passes through the origin of the scat- 
tering potential which means that the three kinks 
collide at one point simultaneously. 

3. The three-particle model explains why the REE 
effect is more pronounced for the solitons collid- 
ing with a small relative velocity. For the particle 
moving in the (a, /3)-plane along the separatrix [red 
line in Fig. [71(c)] , any perturbation results in expo- 
nential in time deviation from the potential ridge. 
High-speed collision results in faster passing of the 
scattering potential and the trajectory of the parti- 
cle cannot be considerably changed. The situation 
is opposite for the slow particle, which corresponds 
to the collision of solitons with a small relative ve- 
locity. 

4. The fractal soliton scattering is explained by the 
presence of the periodic orbits of the particle in the 



11 



scattering potential. Periodic orbits of two types 
were found, each of them is responsible for a par- 
ticular scenario of the particle dynamics, and each 
scenario yields a self-similar pattern for the colli- 
sion outcome as a function of a parameter, such as 
the location of the central effective particle (Sec. 

MM- 

5. Periodic orbits are also responsible for the algebraic 
law p ~ L~ a , where p is the probability to observe 
the three-soliton bound state with the lifetime L 
(see Fig. Eland Sec. MM- 

In the weakly perturbed systems the REE is the dom- 
inant effect. However, if the perturbation is not small, 
the energy exchange effect is mixed with radiation and 
possibly with excitation of internal modes. We thus be- 
lieve that the net effect of inelasticity of soliton collisions 



can be decomposed into three major parts: the radiation- 
less energy exchange, excitation of the soliton's internal 
modes, and emission of radiation. This highlights the 
need for a systematic study as a function of increasing 
deviations from the integrable regime of the relative role 
of these three complementary mechanisms. Such a study 
would be of particular interest for future investigations. 
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